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SUMMARY 


A  siqnif leant  reduction  in  total  test  times  of  computer- 
controlled  fatigue  tests  has  been  achieved  by  the  design  of  a  minimum 
excursion-time  control  system.  The  computer  calculates  in  real  time 
the  minimum- time  path  between  two  excursion  points,  subject  to 
variable  velocity  and  acceleration  constraints.  This  technique  has 
been  extended  to  multi-channel  systems. 
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NOTATION 


Single  Channel  System 
x(t)  time  path 

T  total  excursion  time 

L  total  excursion  distance 

v  maximum  allowable  velocity 

A  maximum  allowable  acceleration 


Multi  Channel  System 

time  path  for  the  ith  channel 
normalised  single  channel  time  path 
total  excursion  time 

total  excursion  distance  for  the  ith  channel 
maximum  allowable  velocity  for  the  ith  channel 
maximum  allowable  acceleration  for  the  ith  channel 


x^t) 

x(t) 

T 


N 

the 

number  of 

channels  with 

non 

L* 

min 

L.  / 

0;  i  -  1,N} 

V* 

min 

{k*  x  v  ; 

i  =  1,N,  L.  f 

0) 

Li  i 

1 

A* 

min 

x  Ai? 

i  =  If N,  Ll  * 

0} 
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1.  INTRODUCTION 

A  minimum  excursion-time  control  system  has  been  implemented 
in  order  to  obtain  a  significant  reduction  in  total  test  time  of  computer- 
controlled  fatigue  tests,  as  well  as  achieving  greater  flexibility  and 
control  over  the  application  of  the  loads. 

Frequently,  load  spectra  used  for  fatigue  tests  are  specified 
as  a  series  of  turning  point  values.  Loads  are  applied  to  specimens  bv 
generating  smooth  waveforms  to  join  successive  turning  points. 

Traditionally,  the  method  of  waveform  generation  has  involved 
the  scalino  of  a  constant  tabular  function,  e.g.  a  sine  function,  which 
has  a  fixed  number  of  steps. 

There  are  several  disadvantaaes  with  this  method: 

(i)  a  constant  tabular  function  leads  to  a  waveform 
which  is  not  optimal  for  excursions  between  most 
turning  points, 

(ii)  a  fixed  number  of  steps  between  turning  points 

causes  load  cycles  of  small  amplitude  to  require 
a  disproportionate  amount  of  time, 

(iii)  there  is  no  independent  control  of  acceleration. 

The  need  for  an  optimal  excursion  aenerator  becomes  apparent 
in  a  situation  where  a  wide  range  of  load  spectra  is  arnlied  to  specimens 
in  tests  of  lonq  duration,  ranging  from  several  weeks  to  several  years. 

The  system  described  below  has  the  following  features, 
which  overcome  the  above  disadvantages: 

(i)  All  calculations  are  performed  in  real-time, 
and  no  stored  function  table  is  required. 

(ii)  Full  control  by  the  operator  of  the  velocity 
and  acceleration  of  each  individual  channel 
exists  (subject  to  physical  limits) . 

(iii)  The  excursion  time  between  any  two  successive 
turning  points  is  always  a  minimum  (subject  to 
velocity  and  acceleration  constraints) . 
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2.  ANALYSIS 


2.1 


manner: 


Single  Channel  System 

The  mathematical  problem  may  be  stated  in  the  following 


What  is  the  minimum-time  path  joining  a  source  and 
destination  such  that 

(i)  the  absolute  velocity  and  absolute  acceleration 
are  bounded, 

(ii)  the  initial  and  terminal  displacements  are 
given, 

(iii)  the  initial  and  terminal  velocities  are  equal 
to  zero. 


The  problem  stated  mathematically  is 


Min  {T} 

such  that  |*(t)|$Vj  |x'(t)j^A;  t  e 
x(o)  «  0;  x(o)  *  0 
x(T)  -  L;  X(T)  «  0  (L>0) 


The  calculus  of  variations  may  be  applied  to  this  problem 
to  prove  that  the  acceleration  may  take  only  one  of  three  values  at 
any  time.  These  values  are  +A,  0,  -A.  Furthermore  it  may  be  shown  that 
the  acceleration  can  only  change  sign  once.  As  this  result  is 
intuitively  clear  the  proof  has  been  omitted. 


There  are  two  possible  solutions  to  this  problem  according 
to  the  relative  sizes  of  L,V,A.  The  first  solution  has  the  acceleration 
taking  on  the  values  +A,  0,  -A  in  sequence.  The  second  solution  has 
the  acceleration  taking  on  the  values  +A,  -A  in  sequence.  Satisfying 
continuity  of  displacement  and  velocity  as  well  as  initial  and  final 
conditions  implies  the  following  solutions. 

1st  Solution  k  >  — 

V  A 


x (t)  -  h  At2 

*(t)  -  At 
X*(t)  -  A 


j 


) 
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x(t)  -  -h  —  +  Vt 
A 

kit)  -  v 
x*  (t)  -  0 


) 

) 

}  V  L 

)  A  V 


x (t)  -  L  -  l*A(t  -  -  -  -) 2 


V  A 


.  L  V 

X(t)  *  -A(t  -  ~  -— ) 
V  A 


x  (t) 


X  ■  —  +  — 
V  A 


L  L  V 

—  <  t  <  ~  +  - 

V  V  A 


2nd  Solution  —  v<  — 

V  A 

x(t)  =  *i  At^  ^ 

) 

x(t)  »  At  ) 

•  •  \ 

x  (t)  -  A  ’ 


x  (t)  -  L  -  *sA(t  -  2 


»<  t  s< 

/  A 


)2  ) 


x(t)  “  -A(t 

x* (t)  *  -A 
T  «  2 


-vf  >  ;  /F*^/F 
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2.2  Multichannel  System 

In  the  multichannel  system  each  hydraulic  actuator  must 
move  in  phase  with  every  other  actuator.  It  follows  that  there  is  an 
equivalent  single  channel  system. 

Since  x  (t)  ■  ~  x (t) ;  x(T)  -  L*,  i  ■  1,N 
1  L 

the  constraints  |x^(t)|  <  i  ■  1,N 


1*1^1  <  Ai 


i  -  1,N 
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imply  )  x  ( t)  j  N<  min  (L*  .  Vi;  i  =  1,N) 

^i 

=  V* 

and  |x(t)  I  <  min  (L*  .  i  =  1,N) 

L. 

1 


=  A* 


L*  could  be  any  reference  length  however  due  to  scaling  considerations 
in  the  computer  implementation  L*  is  defined  as  min  (L^;  i  =  1,N) . 
Substituting  these  values  of  L* ,  V*,  A*  for  L,V,A  in  the  single  channel 
system  determines  x(t),  x, (t)  is  then  determined  from  the  relation 

Xi(t)  =  x(t) 

L* 

3.  COMPUTER  IMPLEMENTATION 

The  system  described  above  has  been  implemented  on  a 
PDP-11  computer  using  the  MACRO- 11  assembly  language. 

The  required  load  spectrum  is  stored  on  a  magnetic  tapd 
in  the  form  of  sets  of  turning  point  values  for  each  actuator,  and  read 
into  memory  in  blocks  as  required. 

The  normalised  single  channel  time  path  x(t)  is  used 
to  generate  a  voltage  excursion  in  the  range  -10  volts  to  +10  volts, 
which  is  output  via  12-bit  digital  to  analoque  converters.  The  signal 
is  ultimately  transmitted  via  servo-hydraulic  actuators  to  produce  the 
required  load  spectrum  on  the  test  specimen,  which  mav  range  from  a 
simple  linear  bar  to  a  complete  airframe. 

A  programmable  real-time  clock  interrupt  is  generated 
every  millisecond,  and  this  is  used  as  the  time  increment  at  which  each 
new  value  of  x(t)  is  calculated.  The  number  of  discrete  values  of  x(t) 
thus  varies  for  each  excursion  according  to  the  values  of  L,  V  and  A. 
However,  during  any  particular  excursion,  the  sane  number  of  discrete 
values  of  x(t)  is  calculated  for  each  actuator  to  ensure  that  thev  all 
reach  their  destinations  simultaneously. 

Due  to  mechanical  and  hydraulic  limitations,  it  has  been 
found  necessary  to  reduce  the  velocity  and  acceleration  of  the  actuators 
if  the  excursion  involves  a  zero-cross inq,  i.e.  aoes  from  tension  into 
compression  or  vice-versa.  Facilities  to  do  this  have  been  incorporated 
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into  the  software  by  allowing  operator  entry  of  maximum  allowable  veloci¬ 
ties  and  accelerations  for  both  zero-crossing  and  non-zero  crossing 
excursions . 


As  floating  point  hardware  is  not  available  in  our  system, 
it  has  been  necessary  to  use  integer  arithmetic  throughout,  to  allow  the 
calculations  to  be  performed  in  real  time.  However  a  hardware  multiply/ 
divide  and  multiple  shift  unit  is  available,  without  which  the  system 
could  not  have  been  implemented. 

For  example,  a  fast  square  root  routine  was  developed, 
based  on  the  well-known  Newton-Raphson  iterative  method,  viz  to  calculate 
X  »  Sq.Rt.(Z),  then  make  a  guess  for  X  and  call  it  Xj,  then  repeatedly 
calculate  a  new  approximation  Xn+1  according  to  the  relation 

*n+l  -  *  <*n  +  f_) 


until  (Xn+i  -  Xn)  is  less  than  some  acceptable  limit  e.  The  technique 
used  in  the  current  implementation  is  to  count  the  number  of  shifts 
required  to  normalize  Z,  i.e.  the  integer  part  of  log  base  2  of  Z,  (this 
takes  only  1  instruction  on  the  PDP-11)  and  use  this  as  an  index  into  a 
table  of  pre-calculated  first  approximations,  i.e.  Xi's  for  the  Newton- 
Raphson  method.  The  correct  (integer)  result  is  then  guaranteed  in  2 
iterations  for  any  16-bit  integer. 

The  algorithm  implemented  in  the  excursion  generator  is 

outlined  below: 

(i)  Input  the  destination,  i.e.  next  turning 
point,  values  for  each  actuator. 

(ii)  Calculate  L*  =  min  {L^;  i  =  1,  n)  where 

Lj_  =»  total  excursion  distance  of  the  ith 
actuator. 

If  the  destination  value  is  equal  to  the  source 
value  for  any  actuator,  then  =  0  and  no 
excursion  is  required  for  that  actuator. 

(iii)  For  each  actuator  determine  whether  or  not  a 
zero-crossing  occurred,  and  set  up  maximum 
allowable  velocities  and  accelerations 
accordingly. 


(iv)  Determine  the  excursion  direction  for  each 
actuator . 
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(v) 

(vi) 


min 

X  V .  ; 

i  =  1,  N; 

Li  / 

0} 

min 

{  —  X  A.; 
L,  1 

i  =  1,  N, 

Li  * 

0} 

(vii)  Calculate  L/V  and  V/A  and  determine  which  condition 
is  satisfied  by  this  case,  i.e.  if  L  N  V  then 

solution  1  is  required  V  A 


if  —  ^  —  then 
V  A 

solution  2  is  required. 


(viii)  A  clock  generating  interrupts  each  millisecond  is 
then  used  to  time  the  excursion  from  source  to 
destination,  with  t  starting  from  0  at  the  start 
of  each  excursion. 


(ix)  The  appropriate  equation  for  x(t),  depending  on  t 
and  whether  solution  1  or  solution  2  is  required, 
as  described  in  the  analysis  section,  is  then 
solved. 


Starting  at  each  millisecond  interrupt,  and  before 
the  next  interrupt  occurs,  the  following  actions 
must  be  carried  out,  until  t  has  reached 

(il  +  Y.)  (case  1)  or  2  /P(case  2): 

V  A  •'A 


(x)  Solve  the  appropriate  equation  for  x(t). 

Li 

(xi)  Multiply  x(t)  by  to  regain  the  correct  scaling 
for  each  actuator. 

(xii)  Convert  the  scaled  x(t)  values  for  each  actuator 
to  units  suitable  for  output  to  the  respective 
digital  to  analogue  converters,  with  appropriate 
magnitude  and  direction  information  included. 

(xiii)  Add  the  base  (source)  value  to  the  above  value 
for  each  actuator. 

(xiv)  Output  the  final  values  to  the  corresponding  D/A 
converters  for  each  actuator. 

Thus  all  actuators  perform  smooth  excursions,  reaching 
their  destinations  together,  independent  of  individual  magnitudes  and 
directions. 
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The  excursion  generator  described  above  requires  approximate¬ 
ly  500  16-bit  words  to  implement,  and  is  incorporated  in  a  large 
general-purpose  fatigue  testing  program  currently  running  on  several 
POP- 11 's  at  ARL. 

4.  CONCLUSION 

A  minimum  excursion-time  control  system  has  been  implemented 
for  the  application  of  load  spectra  to  fatigue  specimens. 

The  system  was  realized  using  a  single-channel  servo- 
hydraulic  actuator  controlled  by  a  PDP-11  computer. 

The  time  taken  to  complete  a  typical  test  was  reduced  from 
6  weeks  to  1  week,  i.e.  a  time  saving  of  approximately  85%. 

While  this  problem  has  been  discussed  in  relation  to  fatique 
testing,  the  method  described  has  wider  applications. 
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